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Abstract 

The Kadomtsev-Petviashvili (KP) nonlinear wave equation is the three dimensional generaliza- 
tion of the Korteveg-de Vries (KdV) equation. We show how to obtain the cylindrical KP (cKP) 
and cartesian KP in relativistic fluid dynamics. The obtained KP equations describe the evolution 
of perturbations in the baryon density in a strongly interacting quark gluon plasma (sQGP) at 
zero temperature. We also show the analytical solitary wave solution for the KP equations in both 
cases. 



1 



I. INTRODUCTION 



The Kadomtsev-Petviashvili equation (KP equation) [l| is a nonlinear wave equation in 
three spatial and one temporal coordinate. It is the generalization of the Korteweg-de Vries 
(KdV) equation to higher dimensions. The KP describes the evolution of long waves of small 
amplitudes with weak dependence on the transverse coordinates. 

The KP equation has been found with the application of the reductive perturbation 
method 2!] to several different problems such as the propagation of solitons in multicompo- 
nent plasmas and dust acoustic waves in hot dust plasmas sl 14|. 



The main goal of this work is to apply the reductive perturbation method |2|-UJ to 



relativistic fluid dynamics 



15 



16| in cylindrical and cartesian coordinates to obtain the 



KP equation. We find that the transverse perturbations in relativistic fluid dynamics may 
generate three dimensional solitary waves. 

Relativistic hydrodynamics is often applied to the study of quark matter in astrophysics 
and heavy ion collisions. Therefore we shall consider an equation of state (EOS) derived 
from QCD jl^]. The obtained energy density and pressure contain derivative terms and a 
wave equation with a dispersive term such as KdV or KP emerges from the formalism. In 



18| . we have performed a similar study in one dimension and found a KdV equation. The 



present work is an extension of |18[ to three dimensions. 

Previous studies on one-dimensional nonlinear waves in cold and warm nuclear matter 
can be found in 19 



This text is organized as follows. In the next section we review the basic formulas of 
relativistic hydrodynamics. In section III we derive the KP equation in detail. In section 
IV we solve these KP equations analytically and in Section V we present some conclusions. 



II. RELATIVISTIC FLUID DYNAMICS 



For a detailed study in relativistic hydrodynamics we su gges t the references |l^ Q. 



The relativistic version of the Euler equation 

dv 



15 



16 



18 



23| is given by: 



dt 



+ {v- V)v 



1 ^dp 



and the relativistic version of the continuity equation for the baryon density is 

dujB" = 



15|: 



(1) 



(2) 



Since Jb'^ = Pb the above equation can be rewritten as |18l. l23l|: 



dpB 2 ( 9v ^ ^ \ ^ , 







(3) 



where 7 = (1 — f ^) is the Lorentz factor. In this work we employ the natural units = 1, 
c= 1. 

Recently we have obtained an EOS for the strongly interacting quark gluon plasma 
(sQGP) at zero temperature. We performed a gluon field separation in "soft" and "hard" 
components, which correspond to low and high momentum components, respectively. In 
this approach the soft gluon fields are replaced by the in-medium gluon condensates. The 
hard gluon fields are treated in a mean field approximation and contribute with derivative 



terms in the equations of motion. Such equations solved 



and space dependence of the quark (or baryon) density 



massless quarks. The energy density is then given by 



17 



properly may provide the time 



26l |. We consider a system with 



/ 27^^ 



/ 27^^ 



and the pressure is: 

/ 27^?2 



V = 
( W 



PB 



PB^'^PB + 

VVBV'(VVB) + i3QCD + 3 



PbV\V'pb) 



Vl6mG^ 



Vps ■ VpB + 



pB ^"^Pb - 



V^PB ^^PB - 



mo 
99' 



27r2 4 

PB V\V^pb) 
V'PB V^(VVb) 



(4) 



■imc 



99' 



ypB ■ y{y^pB) 



-Bqcd + —,— 



(5) 



In (j3]) and ([5]) 7q is the quark degeneracy factor 7q = 2 (spin) x 3 (flavor) = 6 where kp 
is the Fermi momentum defined by the baryon number density: 



PB 



7q 
67r2 



(6) 



The other parameters g, mc and Bqcd are the coupling of the hard gluons, the dynamical 
gluon mass and the bag constant in terms of the gluon condensate, respectively. 
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III. THE KP EQUATION 



We now combine the equations ([T]) and ([3]) to obtain the KP equation which governs 
the space-time evolution of the perturbation in the baryon density using the EOS given 
by (jl]) and As mentioned above, we use the reductive perturbation method |2|-|9| 
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14i\ . Essentially, this formalism consists in expanding both ([T]) and (|3]) in powers of a small 



parameter a. In the following subsections we present the formalism. 

We start with the cylindrical KP (cKP) because it can be converted into an ordinary 
KdV equation (see section V) and similar transverse perturbations considering the two- 
dimensional cylindrical case have also been considered in our previous work {2^. 



A. Three-dimensional cylindrical coordinates 

The field velocity of the relativistic fluid is: 

V = v{r, ip, z, t) = Vr{r, (p, z, t) + ^^(r, ip, z, t) + vl{r, ip, z, t) 



and so \v\ = yfr^ + + v^"^ . 

Rewrite the equations ([T]) and ([3]) in dimensionless variables. The perturbations in baryon 
density occur upon a background of density po (the reference baryon density). It is convenient 
to write the baryon density as the dimensionless quantity: 



p{r, ip, z, t) = (7) 

Po 



and similarly the velocity field as: 



v = ^ (8) 



where is the speed of sound. The components of the velocity are: 

. f ,^ Vr{r,ip,z,t) . v^{r,ip,z,t) 
Vr{r, V?, z, t) = , v^{r, ip, z, t) = 



and 



Vz{r, if, z, t) = (9) 

Co 



We next change variables from the space {r,(p,z,t) to the space (i?, T) using the 
"stretched coordinates" : 



where L is a typical scale of the problem. 

The next step is the expansion of the dimensionless variables in powers of the small 
parameter a: 



,4/3 



1/3 



p = 1 + api + (7 p2 + o- P3 



Vr = CTf n + cr'^Vr^ + a^Vr^ + . . . 
= cr^^^v^, + + a^/^^g + . . . 

v\ = a^l'^v,, + cT^/'t;,, + o'l'^v,, + . . . 

[1 + ((jpi + o'p2 + . . .)]'^' = 1 + \<ypx + \(y''p2 + . . . 

[1 + (api + (7^2 + ■ ■ 0]'^' = 1 + \^Px + \o''pi + • • • 



(11) 

(12) 
(13) 
(14) 
(15) 

(16) 



Finally we neglect terms proportional to a"" for n > 2 and organize the equations as series 
in powers of a, a^/^ and . 

From the Euler equation ([T]) we find for the radial component: 



+ Svr^/Vo'/'c.^ 



dR 



+ 



9pi 
9i? 



) + -n'/'po'" 



dp2 
dR 



r2 



OR 



+ 



+ 



c/ + 

2^2 



V 



dv. 



ri 



dR 



f27gpo\ dpi 2/3^ 4/3 Pi ^Pi 



^27g'po' 



■imc 



2c/ + 4vr2/Vo^/'c,^ 



Pi 



ri 



9i? 



2Vpo^ 



„ 2 I 2/3^ 4/3 „ 2 



9pi /^^^\93pi 



For the angular component: 



a 



3/2^ 



/2Vpo^ 



+ 



+ 37r2/3po-"c 



4/3„ 2 



(17) 



dv, 



and for the component in the z direction: 



1^ 

T 9$ 







3/2, 



I Q^2/3^ 4/3^ 2 

+ 37r ' po ' 



21 



/2Vpo^ 



) + -r^'/'p,'" 



dpi 
dZ 



(18) 



(19) 
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Performing the same calculations for the continuity equation ([3]) we find: 

9 f dvr 



a- 



dvr^ dpi 

~dR ~ m 



a 



OR 



dp2 dpi 
dR dT 



+ Pi- 



dn 



ri dpi 2 y^^n , ^ri 



dVr 



dv 



Cg Vr-^ 



21 



dR ' dR dR ^ T ^ dZ ' T d^ ^ 

In the last four equations each bracket must vanish independently and so {...} = 0. From 
the terms proportional to a we obtain the identity: 

,2 ^ 2x /0'7^2 ^ 2- 



+ 



1 dv. 







(20) 



A 



(21) 



which defines the constant A and from which we obtain the speed of sound for a given 
background density pq. 



27g2 po2 



37r2/3po4/3 



and also 



Pi = Vr^ 

From the terms proportional to cr^/^ using the A constant we find: 



and 



dv^^ _ 1 dpi 
~dR ~ f~d^ 

dvz^ dpi 



(22) 



(23) 



(24) 



(25) 



dR dZ 

Inserting the results ( l2Tl) . ( |23|) . (l24|) and ( !25|) into the terms proportional to in ( fT7|l 
and (l20l) . we find after some algebra, the cylindrical Kadomtsev-Petviashvili (cKP) equation 



d [dpi 
dR\dT 



{2-c/) f27g'po'\{2c,'-l 



2A 



A 



Cs 

6 



Pi 



dpi 
dR 



+ 



9g^Po^ ]&^p 



_^ Pi 
dR^ 2T 



.Smc^L'^A 

From the second identity of fl^Tl) we may write 



1 ^Vi l aVi 



8mG^ 



(26) 



(27) 
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Inserting ( 127|) in the coefficient of the nonhnear term in (126|) the cKP becomes: 



d j dpi 
dR\dT 



3A 



Pi 



dpi 
dR 



d^Pi Pi/ 
dR^ 2T 



1 d^pi Id^pi 







2T2 a<l'2 2 

Returning this cKP equation to the three dimension cyhndrical space yields: 



d_fdp, dp, 13 ,,f^W^ 
dr\ dt ^""^ dr ^ [2^ ' 3A 



CsPl 



dpi 
dr 



9g^ Po^c, 



dx^ 2t 



Pi 



(28) 



(29) 



which is the cKP equation for the second term of the expansion (fTTi) . the small perturbation 
given by pi = a pi. 



B. Three-dimensional cartesian coordinates 

We now write the field velocity of the relativistic fluid as: 

V = v{x, y, t) = v*^{x, y, t) + Vy{x, y, t) + v^{x, y, t) 



and so \v\ = \Jva/ + Vy^ + Vz"^ . 

We follow the same steps described in the cylindrical case to obtain the KP equation: 
1) Rewrite the equations ([1]) and ([3]) in dimensionless variables: 

^(,,y,,,t) = EshyilA (30) 

Po 



V 



(31) 



and so: 



and 



. / v^{x,y,z,t) ^ Vy{x,y,z,t) 
Vx{x, y, z, t) = , Vy{x, y, z, t) = 



Vz{r,ip,z,t) = (32) 

2) Take the equations ([1]) and ([3]) (now in dimensionless variables) from the space (x, y, z, t) 
to the space (X, Y, Z, T) using the "stretched coordinates" : 

X = ^-{x-Cst), Y = —y , Z = —z , T = ^-Cst (33) 
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3) Perform the expansions of the dimensionless variables: 



p = 1 + crpi + a^p2 + a^pz + . . . 


(34) 


Vx = 0-^x1 + ci^Vx^ + cr^Va;3 + . . . 


(35) 


3/2 1 2 1 5/2 1 


(36) 


3/2 1 2 1 5/2 , 


(37) 




(38) 




(39) 


4) Neglect terms proportional to cr" for n > 2 and organize the equations as series in powers 



of a, (T^/^ and cr^. 

After these manipulations the x, y and z components of Euler equation are: 



a< - 



"27g^_p^\ 2 3^2/3 4/3 2 

9X 



dX 



+ 



+ 



^^^)c,^ + 37rVW/^c, 



dX\ 

dvx2 
dX 



or dx , 



,3/2 _ 



o2 ^2 



/27^W\ dpi 2/3 „ 4/3 Pi dpi 
-'f^^fjpi2c/ + AnyW/\ 



U 2 , 2/3, 4/3, 2] ^Pl , /_9^^P0^^ d'^Pl \ . . - 



^^)c.^ + 37r^/W/^c. 



'"'''dX ' V4mG 

r/97^2,„2 



ax 



+ 



/27^X ^2/3 ,/3 



9pi 



dv, 



dX 



^ =0 



and 



,3/2 _ 



ax 



.2/34/3 



api 



(^^)c.^ + 37r^/3po^/^c.^ 



(9w 



ax 



(41) 



(42) 
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For the continuity equation we obtain: 



Again, in the last four equations each bracket must vanish independently. From the terms 
proportional to a we obtain the same A constant as in the cylindrical case given by fl2T|) . 
the same expression for the speed of sound (122|1 and 



{dv^ _ dpi 
dX dX 

^ ^ dvxi dpi 
dX^ dT^^'dX ^ """^ dX 



du 



dv 



pi = Vxi 



(44) 



From the terms proportional to a^/^ we find 



and 



dvy^ _ dpi 
dX ~ dY 

dv^, dpi 



dX dZ 

In ( 14T1) and ( l42l) we have from the terms proportional to a^: 



dv. 



y2 



dv,. 







(45) 



(46) 



(47) 



dX dX 

Inserting the results ( I2TI) . f j44l) . f j45l) . ( l46l) and ( l47l) into the terms proportional to cr^ in (j40 
and (H3ll . we find after some algebra, the Kadomtsev-Petviashvili (KP) equation 0, 14|: 



d \dpi [ (2-c,^) ^ 27^^po^ ^^ (2c,^-l 

dx\dT 



+ 



2A 



_-2/3„ 4/3 



Pi 



dpi 
dX 



Id^p 



1 , Id^Pi 

2dY^ 2 dZ^ 



(48) 



Inserting f l27p in (148 p . the KP with simplified coefficient for the nonlinear term is given by: 



d \dpi 
dX\dT 



3A 



dpi . 
P'dX^ 



9g'po' 



d^pi 



'^'''+l?^ = (49) 



9X3 2dY^ 2 dZ^ 



Rewriting this KP equation back in the three dimensional cartesian space we find: 



dx\dt 



dx 



7r^/3poV3 
3A 



CsPr 



dpi 



+ 



dx V Sitlg'^A 



d^pi 
dx'^ 



which is the KP equation for the small perturbation pi 
expansion ([3l 



Cs d'^pi Cs d'^pi _ 
2 dy^ 2 dz^ 

(50) 

crpi, the second term of the 



C. Some particular cases 



In the one dimensional cartesian relativistic fluid dynamics we have v = v{x, t) and 
Pb = PB{x,t). Repeating all the steps of the last subsection for one dimension, the reductive 
perturbation method reduces to the formalism previously used in [l9l-|25| and we find the 
following particular cases of (!50|) : 

(/) Neglecting the y and z dependence, the fl5Q|) becomes the Korteweg-de Vries equation 
(KdV) similar to the KdV found in 



^ + c — + 
dt dx 



(2 



27(7^0^ - 1) 7r2/3po4/3 



+ 



9^^ Po^c, 



2A 

dx^ 



A 



CsPl 



dpi 
dx 



(51) 



Taking the limit mc — t- oo we obtain from ( 12 ip and 

A 

and flSTj) becomes: 



1 
3 



— + c — + -c " — 
dt ^'^ dx . ,3 '^'^^^ dx 







(52) 



23|, the so called breaking wave equation for pi 



and we recover exactly the result found in 
at zero temperature in the QGP with the MIT equation of state. 

(//) Neglecting the spatial derivatives in (jl]) and ([5]), we obtain flSTl) reduces to: 



dpi , -^Pi , 

+ Cs— h 



dt 



dx 



(2-c,2) f21g^p,\{2cs^-l 



7r2/3po4/3 



Cs 

6 



2A A 

which is also a breaking wave equation for pi with the po, mo and g dependence in its 
coefficients. 



IV. NON-RELATIVISTIC LIMIT 



The non-relativistic version of the continuity equation ([T]) is given by 15|, ll6| : 

dpB 



dt 



+ V ■ (pbv) = 



and for the Euler equation we have 

dv 
di 



15 



16|: 



P 



(54) 



(55) 
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where p is the volumetric density of fluid matter. In this work we study perturbations for 
baryon density in the sQGP fluid, so we deflne the "effective baryon mass in sQGP": 



p = MpB 



(56) 



which will be determined latter. Substituting f l56|) in f l55|) we flnd the non-relativistic version 
for the Euler equation in sQGP: 

1 



dv 
di 



Mp 



B 



Vp 



(57) 



Performing all the calculations described in the last section for the combination of ( 154|) 
and ( 1571) we flnd the cKP equation in non-relativistic hydrodynamics: 



dr\ dt '"^ dr 



CsPl 



dpi 
dr 







dx^ 2t 



and the KP equation in three-dimensional cartesian coordinates: 



d j dp 
dx\ dt 



1 dpi 



dx 



7r2/3pol/3 



3Mc, 



CsPl 



dpi 



dx 



+ 



9^^ Po 



9^ 



Cs d'^pi Cs d^pi 
2 + 2 



(58) 



(59) 



which are the non-relativistic versions of (129|) and (l50l) respectively. During the derivation 
in both cases we flnd from the terms proportional to a in the Euler equation: 

27g^Po 



M 



(60) 



We end this section mentioning that it is possible to obtain ( 158|) and (1591) directly from 
(129|) and (150|1 respectively, performing the two non-relativistic approximations: 



a) 







(61) 



and 



b) 



A = MpoCs 



(62) 



where A is given by ( 12T|) and by ( l60l) . 
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V. THE ANALITYCAL SOLITON-LIKE SOLUTION FOR CKP AND KP EQUA- 
TION 



In this section we present the analytical soliton-like solution to the cKP and KP equation 
given by (^U^ and fISUl) respectively. The KP equation is an integrable system in three di- 
mensions in the same way as the KdV is in one dimension. We introduce a set of coordinates 
that transforms (129!) in an ordinary KdV, which is a solvable equation, and we also present 
the analytical solution of f l50|) . In order to simplify the notation in equations f l29|l and f l50|) 
we define the constants: 



a = 



'1 - c. 



3A 



8m 



and 



Considering the following coordinates j5|, l9|, [11 

4 = r h z 



and 



r = t 



(63) 



(64) 



(65) 



we have: 



d d 

dr 



g3 Q3 

d 



d_ d_ 



Cs<f^d_ 
2 



and 



pi{r,(f,z,t) pi(^,r) 
Using ( 166|) and ( 167|1 in ( 129|) we find the KdV equation in the (^, r) space: 



9r V 2 / 5^ 5^ 
which has the analytical soliton solution given by: 







. 3(m-3c,/2) 2 
pi(^,r) = sech 



a 



{u - 3c,/2) 



4/3 



^ — UT 



(66) 



(67) 



(68) 



(69) 



The exact analytical soliton solution of fl2^ in three cylindrical coordinates is obtained 
substituting fl65|) in fl69|) : 



. . ,^ 3(m - 3c,/2) 2 
pi[r,(p, z,t) = seen < 



a 



{u - 3c,/2) 
4/3 



r + 2; 



t 



(70) 
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where m is a supersonic parameter which satisfies u > ^ and the phase velocity given by 

2 

u + is angle dependent. 

The exact analj^ical soliton solution of the KP equation ( l50l) is given by [gI, Q, 30 |: 



pi(x,y, z,t) = — —sech'^ 

a 



\U-2cs 



x + y + z-Ut\\ (71) 



4/3 , 

where again we have a supersonic parameter U such that U > 2cs . 

VI. CONCLUSIONS 

We have described in detail how to obtain a KP equation in three dimensions in cylin- 
drical and cartesian coordinates in the context of relativistic fluid dynamics of a cold quark 
gluon plasma. To this end, we have used the equation of state derived from QCD in jl?! . 
The resulting nonlinear relativistic wave equations are for small perturbation in the baryon 
density. 

For the cartesian KP the exact soliton solution is a supersonic bump keeping its shape 
without deformation. The cartesian KP contains some particular cases such as KdV and 



the breaking wave equation already encountered in our previous works [18|, |23|. For the 
cylindrical KP (cKP) we also have an exact supersonic soliton solution which deforms slightly 
as time goes on due to the angular dependence in the phase. 

We conclude that relativistic fluid dynamics supports nonlinear solitary waves even with 
the inclusion of transverse perturbations in cylindrical and cartesian geometry. 
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